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In order to identify the Higgs sector realized in nature, the predictions for Higgs boson masses, 
production cross sections and decay widths have to be compared with experimental results. We 
give a brief overview about computer codes for the evaluation of the properties of charged Higgs 
bosons, mostly focusing on the case of the Minimal Supersymmetric Standard Model (MSSM). 
We briefly review the relevance of the various contributions to the charged MSSM Higgs boson 
mass arising at the one-loop level. 
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1. Introdcudon 

Identifying the mechanism of electroweak symmetry breaking will be one of the main goals 
of the LHC. Many possibilities have been studied in the past, of which the most popular ones 
are the Higgs mechanism within the Standard Model (SM) [1], the Two Higgs Doublet Model 
(THDM) [2], and within supersymmetric (SUSY) models [3], which naturally contain at least two 
Higgs doublets. While the SM contains only one physical (neutral) Higgs boson, models with two 
Higgs doublets contain, besides three neutral Higgs bosons, at least two charged Higgs bosons, 
H^. Any evidence for a charged Higgs boson would be an unambiguous sign of an extended Higgs 
sector, i.e. of physics beyond the SM. 

In order to identify the Higgs sector realized in nature, the predictions for couplings to SM 
fermions and bosons as well as self-coupUngs have to be compared with the experimental results. 
To perform this task precise theory predictions are needed for Higgs production cross sections and 
decay widths (or at best for complete processes). Here we give a brief overview about this kind of 
codes concerning the charged Higgs bosons. A recent overview about codes and tools for SUSY 
can be found in Ref. [4]. 

Within the THDM the Higgs potential contains 14 independent parameters, and the mass of 
the charged Higgs boson, M/f±, is usually treated as a free parameter. Consequently, no code exists 
predicting Mjj± . However, tree-level production cross sections and decay width can be evaluated 
with codes like MadGraph [5], Pythia [6], Sherpa [7] or Herwig [8]. We are not aware 
of codes containing higher-order corrections to charged Higgs boson processes. However, this 
could be done in principle with FeynArts/FormCalc [9, 10], for which the THDM model file 
is available. 

Within SUSY many results for the charged Higgs boson beyond the tree-level are available in 
the Minimal Supersymmetric Standard Model (MSSM), see below. The only code that contains 
higher-order corrections for a charged Higgs boson beyond the MSSM is NMSSMTools [11], 
providing masses, production cross sections and decay widths in the NMSSM, i.e. the MSSM with 
an additional Higgs singlet. In the following we will focus on the MSSM. 

2. The charged Higgs in the IVISSIVI 

The Higgs sector of the MSSM contains two Higgs doublets, leading to five physical Higgs 
bosons. At tree-level these are the light and heavy '^,^-even h and H, the '^■^-odd A and the 
charged H^. At lowest order the Higgs sector can be described besides the SM parameters by 
two additional independent parameters, chosen to be the mass of the A boson. Ma (in the case of 
vanishing complex phases), and tanjS = V2/V1, the ratio of the two vacuum expectation values. 
Accordingly, all other masses and couplings can be predicted at tree-level, e.g. the charged Higgs 
boson mass 

ml^=Ml+Ml. (2.1) 

Mz,w denote the masses of the Z and W boson, respectively. This tree-level relation receives 
higher-order corrections, where the loop corrected charged Higgs-boson mass is denoted as Mh± . 

The charged Higgs bosons of the MSSM (or a more general Two Higgs Doublet Model) have 
been searched at LEP and the Tevatron, and will be searched for at the LHC and the ILC. The LEP 
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searches [12, 13], yielded a bound of Mu± > 80 GeV [14]. The Tevatron is especially sensitive at 
low Ma and large tanj3, extending the LEP bounds in this region of parameter space [15]. Also 
at the LHC the prospects for charged Higgs boson searches are best at large tan jS , reaching up 
to Ma ^ 800 GeV [16-18]. At the ILC, if the charged Higgs is in the kinematical reach, a high- 
precision determination of the charged Higgs boson properties will be possible [19,20]. 

Many computer codes exists to evaluate various properties of the (charged) Higgs boson in the 
MSSM: 

• ROE running from a high-energy scale to the electroweak (EW) scale for SUSY parameter is 
provided, for instance, by SoftSusy [21], Spheno [22], Suspect [23] or Isa jet [24]. 
Subsequently they evaluate Mfj± including (some) one-loop corrections. 

• Three codes exist for the calculation of M/f± and the various decay widths, FeynHiggs [25- 
29], CPsuperH [30] and Hdecay [31]. More details on the evaluations can be found below. 

• Calculations for the charged Higgs production at the LHC have been performed in 

Refs. [32, 33] (see also Ref. [34]), for the production in association with a W boson in Ref. [35] 
and for the production of a H^H^ pair in Ref. [36]. However, no dedicated code exists. The 

calculation of Refs. [32, 33] has recently been implemented into Prospino [37] and is also 
included in FeynHiggs [25-29] (including the A/, corrections, see below). More details on 
recent higher-order corrections can be found in Ref. [38]. For the production at the ILC 
we are only aware of one code for the calculation of e+e^ W^H^ [39]. 

• Event generators can produce events with a charged MSSM Higgs boson: Pythia [6], 
Sherpa [7] or Herwig [8]. They contain only relatively few corrections to the Higgs 
boson mass or decay widths. However, they can be linked to the dedicated spectrum and 
decay calculators via the SUSY Les Houches Accord [40]. 

• Indirect constraints on the mass and couplings of the charged MSSM Higgs boson can be set 
via(e.g.) S-physics observables. One dedicated code for this is Superlso [41], another one 
is described in Ref. [42]. Other codes contain some S-physics observables as an additional 
check on the MSSM parameter space. Examples are MicrOMEGAs [43], CPsuperH [30] 
or FeynHiggs [25-29]. LEP and Tevatron bounds on the charged Higgs boson searches 
are currently implemented into HiggsBounds [44]. 

A "comparison" of the dedicated mass and decay width calculators has to be split up into 
the case of the MSSM with real parameters (rMSSM) and the MSSM with complex parameters 
(cMSSM). Hdecay is purely for real parameters. CPsuperH and FeynHiggs can handle also 
parameters with complex phases (for instance A,, Ai„ M^, jx, . . . ). The following corrections to 
Mu± are implemented into the three codes (we do not go into details concerning the attempt to 
capture three-loop corrections): 

(i) Hdecay: Mh± in the rMSSM: 

Ma is input parameter, and the higher-order corrections to Mh± are evaluated. The calculation 
of Ref. [45] (including corrections up to the two-loop level) has been exteded to the complete 
Higgs potential and thus to Mf]±, including the A^, corrections. 
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(ii) CPsuperH:M//± in the cMSSM: 

Mh± is an input parameter, and the higher-order corrections appear as shifts to the neutral 
Higgs boson masses. A large set of one-loop corrections are supplemented by two-loop 
contributions obtained by RGE improvements, see Ref. [46] and references therein. The 
complex Af, corrections (see below) are included. 

(iii) FeynHiggs: Mh± in the rMSSM: 

Ma is an input parameter, and the higher-order corrections to Mfj± are evaluated. They 
comprise the full set of one-loop contributions [28] and various two-loop corrections as sum- 
marized in Ref. [27]. 
FeynHiggs: Mh± in the cMSSM: 

Mh± is an input parameter, and the higher-order corrections appear as shifts to the neutral 
Higgs boson masses. The calculation comprises the full set of one-loop contributions [28] 
and the ^(aftts) corrections as given in Ref. [29]. 

An overview about the various charged MSSM Higgs boson properties implemented into 
FeynHiggs, CPsuperH and Hdecay is shown in Tab. [l]^ The evaluation of the charged Higgs 
decay contains in all three codes the contribution coming from A^ (see below). However, the ef- 
fect of complex phases in the various calculations is taken into account only in FeynHiggs and 
CPsuperH. The Higgs self-couplings are understood to contain an H^H pair. The LHC produc- 
tion cross section calculation in in FeynHiggs is based on Refs. [32, 33], supplemented by the Ah 
corrections. We also list BR(f — > H^b), which be used to evaluate the charged Higgs production 
cross section for Mjj± < nit at the Tevatron or the LHC. 



3. Relevance of higher-order corrections to M^± 



3.1 The higher-order corrections 



As mentioned above, the charged Higgs boson mass, see Eq. ([2.1[), receives higher-order cor- 
rections, which are (at different levels of sophistication) implemented into public codes. In a first 
step in Ref. [48] leading corrections to the sum rule given in Eq. ( p] ) have been evaluated. Then 
one-loop corrections from t/b and t/b loops have been derived in Refs. [49,50]. These were ex- 
tended to the one-loop leading logarithmic terms from all sectors of the MSSM in Ref. [51]. The 
first full one-loop calculation in the Feynman-diagrammatic (FD) approach has been performed in 
Ref. [52], and re-evaluated more recently in Refs. [28,53]. Within the FD approach the leading 
(^(aftt^) two-loop terms have recently been obtained in Ref. [29,54]. In the following we will 
focus on the one-loop corrections. 

Within the FD approach the charged Higgs boson pole mass, M^± , is obtained at the one-loop 
level by solving the equation 

p^-mj,^+t'^X-iP^)=^- (3-1) 



'The Higgs-self couplings and the BR(f — + H'^b) are internally calculated in Hdecay, but are so far not part of the 
output [47], 
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Table 1: Overview about the various charged MSSM Higgs boson properties implemented into 
FeynHiggs, CPsuperH and Hdecay. / denotes a SM fermion, /i, (i = 1,2,3) are the three neutral 
Higgs bosons, denotes a scalar fermion, and xf, xf the neutralinos and charginos, respectively. The 
evaluation of the charged Higgs decay contains in all three codes the contribution coming from Ai,. The effect 
of complex phases in the various calculations is taken into account only in FeynHiggs and CPsuperH. 
Higgs self-couplings are understood to contain an H+H~ pair The LHC production cross section calcula- 
tion in in FeynHiggs is based on Refs. [32, 33], supplemented by the corrections. The BR(? — > H'^b) 
can be used to evaluate the charged Higgs production cross section for M^± < rrit. 



^ff+ff- (^^) denotes the renormalized one-loop charged Higgs boson self-energy. Details about the 
calculation and the full one-loop evaluation can be found in Refs. [28,53]. 

Another class of higher-order corrections that is formally of two-loop order is normally in- 
cluded already into the one-loop result. These corrections originate from the contributions to Higgs 
boson self-energies from the bottom/sbottom sector enhanced by jx tan/3. Large higher-order ef- 
fects can in particular occur in the relation between the bottom-quark mass and the bottom Yukawa 
coupling at large tanjS. Because the tan j3 -enhanced contributions can be numerically relevant, an 
accurate determination of hb from the experimental value of the bottom mass requires a resum- 
mation of such effects to all orders in the perturbative expansion, as described in Refs. [55,56]. 
Effectively the bottom Yukawa couphng is given by 



2Mw l+Af, 



with nib ^ m°^e) = mf\Q) (l + ^) • (3-2) 



nib denotes the running DR bottom quark mass at the mass scale Q, where includes the SIVl 
QCD corrections. The leading one-loop contribution Af, in the Umit of Msusy ^ mt and tan j3 » 1 
takes the simple form [57] 

Af, = ^M3>*tanj3 x/(mg^,m^^,|M3|) + ^A;M*tanj3 x/K,m?,,|M|) , (3.3) 
where the function / arises from the one-loop vertex diagram and scales as 

I{a,b,c) ~ l/max{a^,b^,c^). Here mt^,mf^ and m^^^^nif^ denote the two scalar top and bottom 
masses, respectively. IM3I is the gluino mass parameter, 11 is the Higgs mixing parameter, and 
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At denotes the trilinear Higgs-stop coupling. This type of higher-order corrections being formally 
beyond the one-loop level are included into the Yukawa couplings appearing in the one-loop result. 

3.2 Numerical analysis 

The higher-order corrected Higgs-boson sector has been evaluated with the help of the Fortran 
code FeynHiggs [25-28]. The goal for the theory precision in Mu± should be the anticipated 
experimental resolution or better, where an accuracy of ~ 1.5 GeV at the LHC and ~ 0.5 GeV at 
the LLC could be achievable (see Ref. [54] for a more detailed discussion). 

In Fig. |I| we show AM//± := Mh± — mH± with Mh± evaluated at the one-loop level in various 
steps of approximation. The results are obtained in the no-mixing scenario [17,58]. The upper plot 
in Fig.[l] shows AMh± as a function of tanjS, while the lower plots depicts the variation with jj.. 
Both parameters enter the definition of Ah and are thus potentially important for a precise Mh± 
prediction. In both plots the solid lines represent the full one-loop result including the resum- 
mation, see Eq. The first approximation to this is shown as short-dashed line, where only 

the contributions from SM fermions and their SUSY partners (i.e. all squarks and sleptons) are 
taken into account, still including the corrections. The next step of approximation is shown as 
dot-dashed lines, where only corrections from the t/b and t/b sector are included, still with the Ah 
resummation. The penultimate step of the approximation is to leave out the Ah corrections, but us- 



ing Jfih (i.e. including the SM QCD corrections, see Eq. (3.2)) in the Higgs boson couplings, shown 
as the long-dashed lines. The final step in the approximation is to drop the SM QCD corrections, 
i.e. replacing Mh by nih in the Higgs Yukawa couplings, shown as the dotted lines. 

The dependence on tanjS is analyzed in the upper plot of Fig. |I| We show AMu± for Ma = 
200 GeV and pL = 1000 GeV. The sign and size of the one-loop correction to Mh± depends strongly 
on tanjS, which appears the Higgs couplings to (s)fermions as well as in the Ah corrections. Neg- 
ative corrections are reached for tanj3 ^ 10 with AM//± — 10 GeV for tanj3 0.6. It should be 
kept in mind that within the no-mixing scenario values of tanjS around 1 are excluded by LEP 
Higgs searches [59,60]. Positive values of AMh± are obtained for large tanjS values, reaching 
AMh± « 2 — 6 GeV for tan j8 = 50. The effect of the non-sfermion sector (short-dashed lines) 
is small and stays below 0.5 GeV. The Yukawa coupling independent effects (dot-dashed lines) 
are ~ 2 GeV, largely independent of tan j8 . The contribution from the Ah effects is negligible for 
tan j3 < 5 and grows with increasing tan j8, reaching values of up to 2 GeV for tan j3 = 50. For small 
values of jJ. (not shown here) these corrections stay very small even for the largest tan j8 values. The 
biggest effects on M//± can arise from the inclusion of the SM QCD corrections to nih for tan j3 ^ 5, 
reaching up to 5-10 GeV. 

In order to analyze the dependence of the Mh± prediction on jx we show in lower plot of Fig. |l|. 
AMh± as a function of n for Ma = 200 GeV and tanjS = 50. The pure t /b/t/b corrections (dotted 
line) reach 12-15 GeV. Including the SM QCD corrections (long-dashed) strongly reduced the 
effect to the level of 4-6 GeV. In the next step the A;, effects are included (dot-dashed line). Due to 
Afc oc jU tan j3 the inclusion of Ah results in a strong asymmetry of Mfj± with a large correction for 
negative /i (corresponding to an enhanced bottom Yukawa coupling) and a much smaller correction 
for positive jj. (corresponding to a suppressed bottom Yukawa coupling). AM//± now ranges from 
~ 14 GeV for jx = -1500 GeV to ~ 1 GeV for /i = +1500 GeV. Adding the corrections from 
the other (s)fermions yields an nearly /i -independent upward shift of ~ 2 GeV. Including the 
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Figure 1: AM//± ^//± — is shown as a function of (upper plot), tan j3 (middle) and /i (lower) in 
the m™^" scenario. Mfj± is evaluated at the two-loop (solid) and the one-loop level (dashed). 
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non-(s)fermionic corrections result in another small upward shift of ~ 0.5 GeV. the overall one- 
loop effect ranges from AM//± w 19 GeV for /i = -1500 GeV down to AM//± 3.5 GeV for /i = 
+1500 GeV. 

The results shown in Fig. |l] clearly indicate that for a precise M//± prediction all contributions 
at the one-loop level have to be taken into account. Effects at the GeV level may be probed at the 
LHC and the ILC, see above. Thus, all parts of the one-loop calculation are relevant and should be 
taken into account in a precision analysis of the charged MSSM Higgs boson. 

4. Summary 

In order to identify the Higgs sector realized in nature, the predictions for Higgs boson masses, 
production cross sections and decay widths have to be compared with experimental results. We 
presented a brief overview about computer codes for the evaluation of the properties of charged 
Higgs bosons. While only few codes exist for the THDM or the NMSSM, many different tools 
are available for the MSSM. These comprise tools for RGE running, the evaluation of the charged 
Higgs boson mass, production cross sections and decay widths as well as event generators. Several 
codes provide indirect constraints on the charged Higgs boson sector via the evaluation of B-physics 
observables that can be checked against existing experimental data. We briefly compared the three 
codes available for the evaluation of the mass and decay properties of the charged MSSM Higgs, 
boson: FeynHiggs, CPsuperH and Hdecay. Finally, we briefly reviewed the relevance of the 
various contributions to the charged MSSM Higgs boson mass arising at the one-loop level. It was 
shown all parts of the one-loop calculation are relevant and should be taken into account in any 
precision analysis involving Mu± . 
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